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ABSTRACT 

We investigate two high energy radiation mechanisms, the proton synchrotron and 
the electron inverse Compton emission, and explore their possible signatures in the 
broad-band spectra and in the keV to GeV light curves of gamma-ray burst afterglows. 
We develop a simple analytical approach, allowing also for the effects of photon-photon 
pair production, and explore the conditions under which one or the other of these 
components dominates. We identify three parameter space regions where different 
spectral components dominate: (I) a region where the proton synchrotron and other 
hadron-related emission components dominate, which is small; (II) a region where the 
electron inverse Compton component dominates, which is substantial; (III) a third 
substantial region where electron synchrotron dominates. We discuss the prospects 
and astrophysical implications of directly detecting the inverse Compton and the 
proton high energy components in various bands, in particular in the GeV band with 
future missions such as GLAST, and in the X-ray band with Chandra. We find that 
regime II parameter space is the most favorable regime for high energy emission. The 
inverse Compton component is detectable by GL^S'T within hours for bursts at typical 
cosmological distances, and by Chandra in days if the ambient density is high. 



Subject headings: gamma rays: bursts - radiation mechanisms: non-thermal 



Introduction 



Camma-ray Burst (CRB) afterglows have been detected mainly at longer wavebands, from 
X-rays to radio, extending up to months after the burst triggers (e.g. van Paradijs, Kouveliotou 
& Wijers 2000). These long-lived afterglow emissions are generally interpreted within the fireball 
shock model by the synchrotron emission of external shock-accelerated relativistic electrons. 
The broadband electron synchrotron spectrum (Sari, Piran & Narayan 1998) has proven to be 
a useful paradigm to study afterglow lightcurves in the X-ray, optical, and radio band, and to 
constrain various unknown parameters. On the other hand, there are other high energy emission 
components, whose role in determining the long wavelength afterglows may be secondary or of 
limited duration, but which may, for some time, dominate the high energy. X-ray to GeV spectrum 
of the afterglows. These include the synchrotron self-inverse Compton emission (IC) of the 
electrons (Meszaros, Rees &: Papathanassiou 1994; Waxman 1997; Panaitescu &: Meszaros 1998; 
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Wei & Lu 1998, 2000; Dermer, Bottcher & Chiang 2000a; Dermer, Chiang k Mitman 2000b; 
Panaitescu & Kumar 2000; Sari &; Esin 2001), the proton synchrotron emission (Vietri 1997; 
Bottcher & Dermer 1998; Totani 1998) as well as some other hadron-related emission components 
(Bottcher & Dermer 1998). The emission from such high energy mechanisms in the GeV band 
may have been detected by EGRET in GRB 940217 (Hurley et al. 1994) and probably also in 
some other GRBs, and will be detectable by next generation GeV 7-ray missions such as GLAST. 
In the X-ray band, the electron IC component has been pointed out as potentially detectable at a 
later time (Sari & Esin 2001). 

Previous afterglow snapshot spectral fits to multi-wavelength data up to X-rays using the 
synchrotron emission component (e.g. Wijers & Galama, 1999; Panaitescu & Kumar 2000; van 
Paradijs, Kouveliotou & Wijers, 2000; Freedman & Waxman 2001) have constrained some of the 
fireball shock parameters (e.g. the electron and magnetic field energy equipartition parameters 
eg, e^), which are found to vary over a wide range. Because of the relative looseness of these 
fits, there is a need to tighten the model constraints. A promising way to do this is by extending 
upward the frequency range over which the snapshot spectra are fitted, and to investigate the 
relative importance and additional constraints imposed by the high energy radiation components, 
over a wider region of parameter phase space, which is the subject of this paper. For this purpose, 
we use a simple analytic approach to describe the various spectral components. Following a brief 
general treatment of the particle distribution and the synchrotron spectrum of the electrons (§§), 
we discuss the proton synchrotron spectral component (§|3|) and the electron inverse Compton 
component (§^, and compare their importance relative to the electron synchrotron emission in 
the frequency regime below the electron synchrotron cut-off. Above this cut-off, the relative 
importance between the electron's IC and the hadron-related components is also discussed. In §^ 
we incorporate the high energy absorption due to 7 — 7 pair production in an analytic manner. 
We then present several examples of the broad band spectra of GRB afterglows within different 
parameter regimes in §^, and discuss the detectability of these high energy spectral components 
within various bands, especially in the X-ray band and the GeV band. Finally we summarize our 
findings in §|7[ 



2. Proton and electron distributions and cooling 

We assume that for negligible radiative losses both electrons and protons are shock accelerated 
to a single relativistic power law distribution of index p. Denoting the particle species by the 
subindex "x" (which indicates "e" for electrons and "p" for protons), after a certain time the 
particle distribution becomes a broken power law depending on the relative ordering between 
the injection minimum energy, 7m, a; 5 and the energy at which the particles cool radiatively in an 
expansion timescale, Jcx- Given a maximum particle energy j^^x at which the acceleration time 
equals the energy loss time, the particle distribution is N{'-fx) oc 7"^ for 'ym,x < Ix < lc,x-, and 
N{-ix) oc 7^^"^ for 7c,a: < 7a: < lu,x, if lm,x < 7c,x (slow cooling regime); and is N{-fx) oc 7^^ for 
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lc,x < Ix < lm,x, and A^(7x) oc ^ for -f,n,x < Ix < lu,x, if 7m,z > 7c,x (fast cooling regime). For 
p > 2, the mean particle energy is = [{p — — '^)\lm,x- Hereafter we assume for simplicity 
the same power law index p for both electrons and protons. Assuming that QpU and C,f.n are the 
number densities of the shock-accelerated power-law protons and electrons, respectively, where Cp 
and Ce are the injection fractions and n is the number density of the pre-shock thermal medium, 
one can define the energy portion contained in the power-law distributed electrons and protons 
with respect to the total energy behind the shock to be eg = Ce(?7^e/^7^p)(7e/^), and = Cp(7p/r), 
respectively, where V is the bulk Lorentz factor of the blast wave. Consequently, the minimum 
energy for the power-law distributed electrons and the protons are Q 

The above equations are for p > 2, which for simplicity will be assumed to be valid. For p 
getting closer to 2, the {p — 2)/(p — 1) factor is no longer precise, and when p = 2, it is replaced 
by ln~^(7„^a;/7m,a;)- However, the imprecise of eq.(|^) at p = 2 does not influence the precise of 
the discussions in §^ (e.g. eq.[^) since both correction factors for the proton and the electron 
components are cancelled out. The cooling energy 7c,x is defined by equating the comoving 
adiabatic expansion time t'^^ ~ r/{Tc) ~ Tt (where F is the bulk Lorentz factor of the blastwave at 
the radius r, t ^ F~^(r/c) is the expansion time in the observer frame) to the comoving radiation 
cooling time t'^ = [{t'gy)"^ + (i/c)~^]~^ = t'sy/i^ + ^)) where t'^y and t'lQ are the synchrotron 
and the IC cooling time, and Y = t'^yjt'iQ = (—1 + ^1 + Arjex/ es) /"^ is the Compton factor 
(e.g. Panaitescu & Kumar 2000; Sari &: Esin 2001). Here redshift corrections have been ignored 
for simplicity, r/ = Min[l, {'^m.,xhc,xY~'^] is the fraction of the particle energy that is radiated 
away via both synchrotron and inverse Compton (Sari & Esin 2001), and is the fraction of 
the magnetic energy density with respect to the total energy density behind the shock, so that 
the comoving magnetic field B' = Fc(327rnmpes)^/^. For rjCx/eB *C 1, one has Y ~ ijex/es <C 1 
and t'^ ~ t'gy, and the IC cooling is not important. Alternatively, when rjex/eB ^ 1, IC cooling 
is important, and t'^ ~ {r]ee/eB)~^^'^t'gy (Sari &; Esin 2001). The critical energy above which 
the particle species x cool in an expansion time is 7c, x — (1 ~l~ ^x) ^ {^"^f^xc/^crXfxB^ t), where 
't'sy = (7x"ixC^)/[(4/3)iTT,xC7^(-B'^/87r)] is used, aT,x = (87r/3)(e^/ma;C^)^ is the Thomson cross 
section for particles of mass rux, and cr^.p/crx^e = (w-e/^i^p)^- A caveat is that, by adopting 
Y = (— 1 + \/l + ^rjex/ es) /2, one has implicitly assumed Ljc/Lgy ~ Usy/Us, where Ugy, Ub are 
the synchrotron radiation and magnetic field energy densities, respectively. Strictly speaking, 
this is valid only when the IC cooling occurs in the classical Thomson regime. In practice, high 
energy electrons may cool in the Klein-Nishina regime, especially at the early afterglow phase. 
The cooling frequency in the above analytic treatment is then a crude approximation if the energy 
where the electron cool is in the K-N regime, and more careful treatment ought to be made. As 



^We note that the larger minimum proton energy estimates 7m,p = T adopted, e.g., by Vietri 1997; Bottcher & 
Dermer 1998 and Totani 1998, when p is close to 2 as assumed in those references, lead to larger proton synchrotron 
flux levels than those obtained using equation (hi), because of the higher 7^ involved. 



- 4 - 



discussed more in §^ the analytic treatment presented in this paper is adequate at one hour or 
even earher after the burst trigger, but may be crude for the prompt phase. 

The maximum energy ^u,x is defined by equating the comoving acceleration time 
= 2TTari^/c = 2TTa'yu,x'mxc/ eB' (where tl is the Larmor radius, and a is a factor of the 
order of unity for relativistic shocks) to the minimum of the comoving adiabatic cooling time 
t'^^ and the comoving radiation cooling time t'^. This gives ^u,x = M.\n.{{2T:a)~^{eB' /mxc){rt), 
[3e/aB'aT,xi^ + ^x)]^'^^}, where for protons the first part of the equation applies, while for the 
electrons it is the latter. 

The electron synchrotron emission spectrum is a broken power law separated by three 
characteristic frequencies: the self-absorption frequency i^a,e, the characteristic frequency for the 
minimum energy particles, Um,e ^ (|r) |f^7m,e) and the coohng frequency, z/c,e ^ Jp|f^7c,e- 
Similar expressions apply for protons, except that synchrotron absorption and cooling are less 
important than in electrons. For both p and e the maximum particle energy 'ju,x defines a 
synchrotron cut-off frequency u^^x — i^^)^'^rc^u,x- the above expressions, the factor (4/3) 
gives a more precise conversion from the comoving frame to the observer frame (Wijers & Galama 
1999). The spectral indices of the four segments ordering from low to high in frequency are [2, 1/3, 
-1/2, -p/2] for the fast-cooling {v^e < ^m,e) regime, or [2, 1/3, -(p-l)/2, -p/2] for the slow-cooling 
(fc,e > i'm,e) regime, respectively (Sari et al. 1998). 



3. Synchrotron components 

A potentially interesting high energy emission component of GRB afterglows is the synchrotron 
emission from the shock-accelerated relativistic protons (Vietri 1997; Bottcher &: Dermer 1998; 
Totani 1998, 2000). Compared to electrons, the protons are inefficient emitters due to their much 
larger mass. The ratio of the characteristic synchrotron frequencies of protons and electrons is 

Vm,p _ ( Tm^y ( _ ( gp/Cp ^ / 'me_ 

) \(^e/CeJ [rrip 

where eq.(|l|) is used. The peak synchrotron fiux for the emission of particle x is 
Fu,max,x oc n^jP^, oc n.j./m^., where P^, = (j)o{V^e^B' /m^c^), and (po ~ 0.6 (Wijers & 

m,x m,x 

Galama, 1999). This is independent of whether the peak is at i'rn,x or at Uc^x- Thus the ratio of 
the peak fiux of the two particle species is 

Pu,max,p f ( 

Fu,max,e \Ce/ \ 

We see that both u^.p and F^-max^p are much smaller than Vrn,e and Fy^rnax,e, respectively, 
indicating that the proton emission component is usually buried under the electron emission 
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component. However, the ratio of the coohng frequencies is 

2 

(4) 



(see which means that protons barely cool while electrons cool rapidly. This allows, in some 
cases, the proton component to dominate over the electron component at high frequencies. 

We consider now a certain observation band Azv around u which is above both the cooling 
frequency and the characteristic synchrotron frequency of the electrons, but is below the electrons' 
synchrotron cut-off frequency, i.e. fc,e(and Vrn,e) < v < i'u,e- Regardless of the relative ordering 
of Vm,e aiid fc,e (slow-cooling or fast cooling), the electron synchrotron flux at this frequency 
is F^^e{i^) = {i^ /^m,e)~^^~^^^'^ F^,max,e ' i'^c,e/^y^^ ■ If < i^u,p is also Satisfied, the flux of the 
proton synchrotron component at the same frequency is Fu^p{v) = {v / Vm,p)~''^^~^^^'^ Fy^rnax,pi since 
Vm,p < Vm,e and usually i^c,p>i^u,p- Thus 

Fu,e{l^) \F^ ,max,e j \^m,e J V ^c,e 

(5) 

We see that v must be ^ u^e to make the proton component show up in the spectrum, and that a 
large Cp and a small Ce favor the proton component. The dependence on the injection parameters 
Ce and (p is weak for p close to 2. 

For the purposes of numerical examples we use below a standard scenario (e.g. Meszaros &: 
Rees 1997) in which a blast wave with total energy per solid angle £ = E/Q = lO^^ergs/Sr £^2 
expands into a uniform interstellar medium of particle number density n cm~^, and J7 is a 
putative jet opening solid angle. (In follows we will not discuss the jet dynamics. This does 
not influence our discussions on the early GeV afterglows, but may quantitatively, although not 
qualitatively, change the discussions of the late X-ray afterglows, if the transitions occur after the 
jet breaks. A wind-like external medium with n oc may be also incorporated in a similar way 
but we do not discuss this here). For the uniform external medium case, the blast wave starts to 
decelerate at a radius r^ec ~ (3<S/47rnmpC^rg)^/^ ~ 2.6 x lO^^cm £l2'^n~^^^TQ^QQ, at an observer 
time tdec ~ (rdec/^rlc) • (1 + z) ~ 2.4s £U^n-y^T-%l{l + z) s. Here Fq = 300Fo,3oo = E/M^c^ 
is the initial bulk Lorentz factor of the shocked materials at t = t^ec and the factor {\ + z) 
reflects the cosmological time dilation effect. After collecting enough interstellar medium 
materials, the blast wave will be decelerated self-similarly and its advance is measured by 
r{t) ~ [12.5ct/(l + z)/47rnmpc2]V4 = i.g x lO^^cm (.552/71)^4 + z)^/^^ while the bulk Lorentz 
factor decays as F(t) ~ [3£:(1 + zf/A^ ■ A-irnmpC^t^]^/^ = 19.4(<S52/n)i/^[t/j/(l + z)]-^/^. Here 
r = 4F^ct(l + z) has been adopted, and th is the earth observer time in unit of hours. The 
comoving magnetic field is B' = 7.5G e]/^i552^n^/^[t/i/(l + z)]^'^/^ G. The breaks in the electron 
energy distribution are then given by 

7„^,e = 5.9 X 103(e,/Ce)(^52/n)V8[t J(l + z)]-3/8, (6) 
7e,e = 1.96 X 102(1 + Y,)-h^'£^2^'n-^/^[th/{l + z)]V8. (7) 
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In eq.(P) and hereafter p = 2.2 for both electrons and protons has been adopted to 
calculate numerically the coefficients. The p-dependences in the power indices, if necessary, 
will be still retained. The observed characteristic electron synchrotron frequency is 
then i^rn,e = 2.9 X 10i6Hz(ee/Ce)^eB^'?52^i/^^^^(l + z)^/'^, while that for the protons is 
t'm.p = 4.6 X 10^Hz(ep/Cp)^e]/^<f52^t/i^^^(l + zY^"^. The electron cooling frequency is 

i^ce = 3.1 X 1013Hz(1 + Y,)-h-//^£^2'^\-\'/\l + (8) 

and the cut-off maximum frequencies for the electrons and the protons are 

u^^e = 2.3 X 1023Hz a-i(l + Y,)-\£,2/ny/%'^\l + z)-'/^ (9) 
u^^p = 2.8 X 1023Hz a-h^ji\£52nf/%^^\l + z)-^/\ (10) 

Finally, the maximum electron synchrotron flux is Fornax e = [(4vr/3)r^(^e'^r-P,y /47rD^l(l + z) = 

i(l)oV3/3){e^Cen/mc^D'^)TB'r^{l + z) = 29mJyCee]i'^n^/'^£52D2i{l + z), where 
D = D{z) = 1028cmD28 is the proper distance of the source (which, depending on the 
cosmological model, is also redshift dependent). Similarly, the maximum proton synchrotron flux 
is F^^rnax,p = 15.7/iJyCpe^^^n^/2^g2Z?^g2(l + z). 

For V < iyu,e, the condition for the proton synchrotron component to overcome the electron 
component is Fy^p{vu^p) > F^^e{'^u,p)- Using equations (^), (§) and (p^), this translates into 

(1 + Y^f/^B > 594(e,/6,)2{*'-l)/3(C,/Ce)^(^-^)/^a2/3£:-//^2^-Vl2[,^/(l + ,)]-l/12^ (11) 

which is shown as line 1 in Fig. 1. 

Another condition for the competition between the p and e components which might be 
thought to be relevant is Uu^p > Vu,e, which corresponds to (1 + Ye)'^/^eB > 0.88a'^^'^£^2^^'^ 
n~^/^-^[t/j/(l + z)]~^/^^. However, this condition may not be essential, since as shown by Bottcher 
&: Dermer (1998) there can be other proton-induced electromagnetic signals extending above this 
energy, of level comparable to or lower than the proton synchrotron component. These include the 
synchrotron radiation from the positrons produced by vr^ decay and the 7-rays produced directly 
from TT^ decay. These components may be regarded as an extension to the proton synchrotron 
component, which would stick out above the electron synchrotron component even if i^^^p < i'u,e- 
However, these proton components will compete with the electron IC component, which we discuss 
in §|. 

4. Inverse Compton component 

The inverse Compton (IC) is mainly important for the electron component, both for the 
electron cooling and for forming a separate high energy emission component^. The condition for 



^The IC of the protons is not important since r/p — {'ym,p/'yc,pY ^ <C 1. 
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the IC cooling to be important is rjEe/eB > 1, which has been exphcitly addressed by Sari &; Esin 
(2001). Here we investigate the condition that the IC emission component overtakes the electron 
and/or the proton synchrotron component. 

The IC component has a similar spectral shape to the synchrotron component, but its low 
energy peak is Lorcntz-boostcd by roughly a factor of 7^6' while the frequency spread where 
it is important is stretched out, by comparison to the spread of the synchrotron spectrum, 
extending between the boosted characteristic and cooling frequencies. Sari Sz Esin (2001) have 
explicitly presented analytic expressions for the IC spectral component, and found that the 
power law approximation is no longer accurate at > since electrons with a range of 

Lorentz factors between j^n and jc contribute equally to the emission at each frequency. For the 
convenience of the following discussions, we will still adopt the broken power law approximation 
to perform order-of-magnitude estimations, bearing in mind that more accurate expressions would 
be necessary in more detailed calculations. 

In this approximation, the IC spectral component can be represented (Sari & Esin, 2001) by a 
four-segment broken power law with power indices ordered from low to high frequency of [1, 1/3, 
-(p-l)/2, -p/2] in the slow-cooling regime, or of [1, 1/3, -1/2, -p/2] in the fast-cooling regime. The 
break frequencies are = Jm,e^a,e, i'm,e — lm,e^m,e, and the IC cooling frequency ujfl ~ lc,e^c,e- 
The maximum flux of the inverse Compton component is roughly a factor {u'pfju'^){i^m.e/i^m,e) 
that of the synchrotron component, where u'^f^ ~ (4/3)c(j7-^e^^7e " (r/Tc) ■ [ATQn) and u'^ = S'^/Svr 
are the comoving synchrotron photon and magnetic field energy densities, respectively. This gives 

- ^^T,eCenr = 3.5 X 10-'Cenn7 = 5.7 x 10-^Ce4^'n3/4[t,/(l + ^)]V4 . (12) 

This shows that generally the IC component can only overtake the synchrotron component beyond 
the synchrotron component's cooling break, but before the IC component's cooling break. In the 
slow-cooling regime, for a frequency u satisfying z^ce < < i'u,e and i/^g < z^ < ^c^, the flux ratio 
of the IC and the synchrotron components is 



F'^(u) /i^IC ^ ^ IC ^(p-l)/2 / xl/2 



Alternatively, in the fast-cooling regime, for a frequency u satisfying e < < ^u,e and 
v^^e ^ ^ ^ ^m,ei ^^^^ ratio of the IC and the synchrotron components is 

FlCf^\ /Z7IC \ /.,IC\l/2/ \ (p-l)/2 



(13) 




(14) 



Similarly to the proton synchrotron case, the conditions that the IC component overcomes the 
synchrotron component is FI^^{v^^) > F^^g{i'l^) for the slow-cooling case, or Fl'^{iy^^^) > -^iy,e(2^^e) 
for the fast-cooling case, and both conditions can be simplified to 

16 

yC7T,eCerir7e,e7fc^^ > 1- (15) 
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Using eqs.(P) and (0), this IC-dominance condition over electron synchrotron can be re-expressed 

(1 + Y,)eB < 3.8(ee/Ce)(^-'^Ce(^52/n)(^-')/^[t„/(l + (16) 



as 



This is the line labeled 2 in Fig. 1. We see that a large eg or a small es makes the IC component 
more prominent, since a large eg enhances 7m, e and a small €b tends to increase 7c,e due to the 
inefficient synchrotron cooling. A denser medium also favors the IC component (eq.||T^). One 
comment is that if one takes into account the flux increase above due to the logarithmic term 
from the scattering contributions of the different electrons (Sari & Esin 2001), the IC dominance 
condition is less stringent than (|lB). 



The cut-off energy in the IC component is defined by = Min(7^ gZ^„ ei ^^^n e)' 
where ^ is the Klein-Nishina limit. A rough estimate of this frequency is given by 
^KN,e ~ 7KN,e^KN.e> ^here 1^^^ ^ ~ it^) f^^^L ,e ^ ^ud 7^^^ ~ ^/Min (/ii/^.^^ , /iz^^ J = 
Max(7KM,,i,7KN,e,2), where 7KN,e,i = 1-1 x lOHCe/eefe'^'/'s-.'^^-y^Til + zr^^ and 
7KN,e,2 = 1-0 X 10^(1 +Ye)'^e^J'^ £^2^ 717^^11/^(1 + zy/^ . In any case, the condition » Vu,e always 
holds, which means that above z^^^e the IC component always sticks out. By comparing both 
7j<,j^ ^ and 7KNe2 with ^c,e (eq.0]) and noticing their temporal dependences, we can see that at 
t/t = 1, both 7j^j^ ^ ^ and 7j^j^ ^ ^ are far greater than 7c,e, and only 7j^j^ ^ ^ may becomes smaller than 
7c,e at the prompt phase. This means that throughout the afterglow phase, the scatterings with 
the peak-flux background photons {v'^ g for slow-cooling and v'^ g for fast-cooling) for the cooling 
energy electrons are well in the classical regime, and at t/i ~ 1 or even earlier, the scatterings with 
the peak-energy- flux background photons [v'^ g for slow-cooling and v'^ g for fast-cooling) for the 
cooling energy electrons are also in the classical regime. This justifies the approximate analytic 
cooling treatment adopted in this paper, at least for th > 1- For earlier afterglows, especially at 
the prompt phase, the cooling frequency estimate may be still viable, but the spectrum above the 
cooling frequency ought to be treated more carefully by taking into account the Klein-Nishina 
modifications. 

In the ee,eB space (Fig.l), we see that the regions defined by (11) and ( |T6| ) do not overlap. 
This means that for frequencies v < Uu^e-, given a certain set of parameters, only one or the 
other high energy component (IC or proton synchrotron) competes with the electron synchrotron 
component. We note that there is also a substantial region of parameter space in which neither 
of these high energy components can be dominant. For frequencies v > Vu^e-, on the other hand, 
there is an overlap of the proton dominated phase space region with the IC dominated phase 
space region in the €e,€B plane. We thus need to compare the relative importance of these two 
components. Below, we choose the slow-cooling case as an example. The fast-cooling case can be 
discussed in a similar way and results in qualitatively similar conclusions. 

Within the slow-cooling regime, there are two cases. If Uu^e < < t'g^, the IC emission 
is important, and one can use eqs. (|5|), (p^) and (13) to define the condition that the proton 
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component overcomes the IC component^, i.e., F,y^p{u) / Fl'^{i') = (ep/ee)^^ ^^(Ce/Cp)^^ 
(me/mp)(3p-i)/2/[(16/3)aT,eCenr7^;'^] > 1. This can be translated into 

. 1 - ^ P±l_ 7-p 3p-5 

ee < 8.1 X 10-^e|Cr'Cp £52^''"''^ n'm^) [t^/{i + 2)] lelF^ . (17) 

This is hne 3 in Fig. 1, to the left of which proton synchrotron dominates over IC, which again is 
a relatively small phase space region in terms of eg (but larger than the region left of line 1, where 
proton synchrotron dominates over electron synchrotron). 

The second case in the slow-cooling regime is for v > f^^, and here the competition 
between the two components is more difficult to quantify. The fact that usually Vu,p i^^e 
would seem to indicate that the proton component is not important. However, as Bottcher 
& Dermer (1998) have shown, photo-meson interactions between the relativistic protons and 
the low energy photon spectrum lead to additional hadron-related spectral components at 
high energies, with a i^Fi, level which can be a substantial fraction of that of the proton 
synchrotron component, but extending to much higher energies. For our purposes it is 
not necessary here to compute these components in detail, assuming instead as a rough 
estimate that beyond the maximum proton synchrotron frequency i'u,p the extended hadron 
component uFi, flux level represents a fraction k = O.lk-i of the proton synchrotron flux level 
at the proton synchrotron cut-off, F^^p{uu^p) = F^^rnax,p{i'u,p/i^m,p)^^^^^^^'^ , so that one has 
F^^p{u)/F^,^ 

.p{i^u,p) — 0.1/c_i (j^u^p/j^). This is consistent with Bottcher & Dermer's Fig.l, and 
in general k = O.l/c-i represents an overestimate of the flux level from these hadron-related 
components. We then compare F,^p{u) to FI%{u) = Fl%,,^,{ul';i/u}^J-(-P-^y^{u/ul'^,)-P/\ 
Again with the help of §), (||) and (||), and noticing FI^^{u) / Fl^^i^l^^) = (u/i'l^e)'^^'^ 
and F^^e{i'u,p)/Fu,e{i^l1) = (^^«,p/^^c*^)"^^^ one can get F^^p{v) / Fl^^{v) = [F^^p{v) / F^^p{vu,p)] ■ 
[F.,p{v^,p)/F,Av^,p)] . [F,A^^,p)/F,,M^)] . [F,A^l%)/Fl^eii^l%)] ' [Fl^eii^l^) / Fl^ei^^)] = {ep/eef^-'^ 
(C/Cp)(^-')("^e/mp)(3p-i)/2/[(16/3)aT,eCenr7^;'^] • {i^u,p/ i^]^eY'\i^ / i^uj^-^^'^ ■ After some 
further derivations, the condition F^p{v) > Fl^^{v) can be translated into 

(1 + Y,f-^'>eB > 215A:Zr'el-02e;0-51<^-0A7^0.09^0.34_j-^0.22^-0.35^-0.10(i + ^)0.01^-0.04^ (18) 

These are the lines labeled 4 (calculated for k = 0.1) and 4' (for A; = 1) in Fig. 1. In the 
above expression we have adopted p = 2.2 for both electrons and protons in order to avoid an 
unnecessarily complicated expression. (The explicit spectral index p-dependence can be written 
out straightforwardly, and is presented in the Appendix). Notice that this criterion is mildly 
dependent on u, and we have adopted v = lO^^Hz z^26 for the typical frequency in ([l^). Above 
this frequency (several TeV), the self-absorption due to 7 — 7 pair production becomes important 



^Here we have assumed a — (p— l)/2 slope for the proton component extending to infinity. When u > Vu.p, this is 
an over-estimate for the other hadron-related components (Bottcher & Dermer 1998). Thus the real criterion is even 
more stringent than (E^. 
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(§1) and the comparison is no longer meaningful. For lower frequencies, the constraint on the 
hadron-dominant region of parameter space is more stringent. 

In summary, the competition between various components in the ee,e_B diagram can be read 
off from Fig. 1. A) The first case, applicable to frequencies v < i'u,e, is shown by the solid lines 
(eqs.|jlT| and |jl^), which divide the space into three regions, in which the proton synchrotron 
component competes with the electron synchrotron component, and the electron synchrotron 
competes with the the electron IC component, respectively. Regions I and II are the regions where 
the proton component and the electron IC component dominate in the spectrum, respectively, 
while in region III, neither of these high energy spectral components can overcome the electron 
synchrotron component. B) The second case, applicable to frequencies v > i'u,e, is shown by the 
dashed lines dividing the parameter space into two regions, in which the proton synchrotron and 
the electron IC components compete with each other. Region I' is where the proton component 
as well as other hadron-related components may overcome the electron IC emission component, 
while region IF is the IC-dominated region. The dashed line 3 (eq.||l^) is the separation line 
for the case of < f^*^, and the dashed line 4 (eq.||l^) is the separation line for the case of 
u > v'l'^, calculated for a typical frequency u ~ 10^^ Hz. For lower frequencies, this line moves 
leftwards, causing the hadron-dominant phase space region to shrink. It can be seen that the 
hadron-related components are usually masked by the electron IC component, unless €e is very 
small. We note that both eqs. (0) and (ll) are derived under conditions which are maximally 
favorable for the proton components, which includes also adopting the extreme case of Cp = 1, 
which may be an over estimate (Vietri 1997). Bearing all these facts in mind, we expect that the 
actual proton-dominated regimes (I and I') could be even smaller than what is indicated in Fig.l. 
On the other hand, in both the low energy band {v < i'u,e) and in the high energy band {v > i'u,e), 
the IC component is important in a much larger portion of the parameter phase space (regions II 
and IF). 



5. 77 pair attenuation 

For the higher energy GRB photons in the observer frame, the comoving frame photon energy 
exceeds mc^, and for sufficient high photon densities, a photon with energy E in the observer 
frame may be attenuated by pair production through interaction with softer photons whose energy 
(also in the observer frame) is equal to or greater than Ean = {Tmc^)'^ /E{1 + z)^, depending on 
the impact angle between the two photons. This may greatly degrade the high-energy fluence 
level, and the corrections due to this 7 — 7 absorption process needs to be taken into account. The 
7 — 7 absorption in the GRB prompt phase has been studied by several authors, e.g. Krolik & 
Pier 1991; Fenimore et al. 1993; Woods & Loeb 1995; Baring & Harding 1997; Lithwick &: Sari 
2001. Here, however, instead of the prompt phase we concentrate on the afterglow phase. 

To treat the absorption in the afterglow phase we adopt an analytical approach similar to the 
one developed for internal shocks by Lithwick & Sari (2001, their eq.[2]), which we adapt here to 



-li- 



the external shock scenario. Instead of using 5T (the temporal variation timescale in the internal 
shock scenario), we use in our case an emission timescale t/(l + z), which is the expansion time as 
viewed by the earth observer with the cosmological time dilation effect correction. Assuming that 
the emission spectrum around Ean = hvan is Ly{v) = Ly^Q{v /van)~^ , the total photon number with 
E > Ean can be estimated as iV>£;„„ ~ j^^{Li,fi/hi'){v /van)~^dv-t/{l + z) = L^{i'an)[t / (1+ z)]/ hp . 
We adopt an averaged 77 cross section given by CaT,e, where C is a constant dependent on the 
photon spectral index. Svensson (1987) gives an analytic expression leading to C = 11/180 for a 
photon energy density index of -1 (photon number density index -2), also used by Lithwick &: Sari 
(2001). More detailed calculations by, e.g., Coppi &: Blandford (1990) and Bottcher & Schlickeiser 
(1997) lead to a slightly larger C of the order 0.1. Noticing that Ly{van) = Fy{yan) ■ 47r 1)^/(1 + z), 
the attenuation optical depth is 

^77l^J 47r[4r2ct/(l + z)]2 lQWhl3t ' ^ > 

where in terms of the quantities t and Fy{van) as measured by an Earth observer, the redshift 
factor has cancelled out. 

It is seen from eq. ( [T^) that the dominant spectral dependence of T-y^(i^) is on Fy{van)- For 
~ 1 TeV (where absorption becomes important), we find Fan ~ 2.6keV T^l + 2)"^, which is 
above /ifc,e- In this band, the electron synchrotron emission component dominates in a large region 
of phase space (Fig.l), and even if the IC component may potentially dominate, this happens 
at a later time (§^, when the GeV-TeV emission is not important. Thus for the regime we are 
interested in, we can approximate Fi,{van) = Fi,^max,e{^c,e/ J^m,,e)~^^~'^^^'^ {van/ Vc^e)"^^"^ , and take 
(3 = p/2. For the model of the dynamics adopted used here and for p = 2.2, eq.(|^) reads then 
T^^{v) ~ 0.56C_i(l + Y,)-\l + z)(7^'-8)/8eg-2)/4eP-ic-(P-2)4+^)/^{P+4)/8i(8-3p)/8^p/2^ ^20) 

where C_i = C/0.1, and the dependence of (1 + l^e)"^ may be dropped for the regime I bursts 
(ee ^ es and ~ 0). Notice the mild dependence on v and the weak dependence on t. The 
absorption becomes important only when v approaches 1 TeV (~ 2.4 x 10^6 Hz). This simple 
treatment is in qualitative agreement with the more detailed simulations of Dermer et al. 2000b 
(their Fig.3). 

An approximate expression for the final spectral flux including 77 attenuation is given by the 
flux escaping from a skin-depth of unit optical depth, or 

Ff{v) = Fi"\u)/{l + T,,\ (21) 

where Fy{vy°^ includes the contributions from all the spectral components discussed above. 



6. Results and implications for high energy observations 

With a simple numerical code that includes the three spectral components discussed in this 
paper, as well as the 77 attenuation effect modeled through eq.([l9|), we have investigated the 
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spectra and the lightcurves in various energy bands for different choices of the most relevant 
parameters, especiahy and es- Because of the relatively small value of r^^ for most energies of 
interest here, we have treated this as a simple absorption process and have neglected for simplicity 
the effects of the secondary pairs which it produces. The results generally confirm the division 
of the €e,€B phase space sketched in Fig.l. Below we present some examples and explore the 
detectability of the proton component and the IC component in the afterglow phase in various 
bands, in particular at GeV energies with future missions such as GLAST, and in the X-ray band 
with the Chandra X-ray Observatory. The signatures of the IC and proton components may be 
detected at high energies in at least two ways. One is through snapshot spectral fits, which require 
a wide and well-sampled energy coverage, including the MeV to GeV band. Such simultaneous 
measurements may be achieved in the Swift and GLAST era. Another, simpler way, is to study 
the lightcurves at some fixed high energy band, looking for a possible hardening of the lightcurve, 
the details of which we discuss below. 



6.1. Detectability of the proton synchrotron component 

To explore the detectability of the proton synchrotron component, we choose as typical 
parameters for the regimes I and V the values es = O.Ses^.s, eg = 10~'^ee,-3 and n = 10^n2. Notice 
that es = 0-5eB,.5, Cg = 10~^ee,-3 would fall outside region I for the density n ~ 1 assumed in our 
Fig. 1, but adopting here a higher density n2 ~ 1 these parameters are appropriate for region I. 
A denser medium will reduce fc,e (eq.[§) and consequently enlarges the region I (eq.0|). We do 
not explore an even lower eg, since it is unlikely to have values ee<rne/mp 0.5 x 10^^. 

For the frequency range ly < Uu,e, the critical time tp at which the proton component overtakes 
the electron component, tp, can be derived by making eq.@ greater than unity. From eq.(|8|), and 
taking a certain band, e.g., u = 10^^ Hz 1^23, one gets 

tp = 1.1 hr e-^%iee,^s/ep)^^^-'\Cp/Ce)'^'''^£,2Wa + ^)"'^23'- (22) 

The dependence on n is steep, so that for a lower density medium the overtaking time tp could be 
too late for observational purposes. The dependence on the frequency is also very steep, so that 
the overtaking time for slightly higher v shifts to much shorter time. For example, for 1^23 = 10 (4 
GeV), the overtaking time moves down to tp ~ 40 s. 



Equation (22) indicates that for the above set of parameters, the proton synchrotron 



component will show up in the GeV band. The possibility of directly detecting this component is 
interesting for several reasons. One is in possibly providing a constraint on the ratio of the proton 
to the electron injection fractions into the acceleration process, Cp/Ce, which is for the interest 
of the shock physics. This would also have important implications for the overall energetics of 
the fireball. A detection might also provide information on whether the index p is the same 
for electrons and protons, as assumed here for simplicity (but departures from which would be 
interesting). We note, however, that the absolute luminosity of this proton GeV emission is too 
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faint to be detected by GLAST, even for bursts at a distance of z ~ 0.1, unless <f52>150. It 
may, however, be detectable for <?52>1 by future larger larger effective area ground-based GeV 
telescopes from some close regime I bursts (see §6.4 and Fig. 4). 

In the X-ray band where Chandra is sensitive (z^ ~ i^is)) because of the steep frequency 
dependence tp oc z^~^ the overtaking time occurs always too late, when the flux is low. Thus it 
appears impossible to detect the proton synchrotron component in the X-ray band. Figure 2a 
shows the snapshot spectra including the X-ray and GeV region for the typical parameter set 
representative of regime I, with a relatively dense external medium {n2 ~ 1). 

One remark is that the condition for the proton synchrotron component to dominate a certain 
high energy band is that fc,e must decrease with time, to make the electron component relatively 
less prominent (see eq. |^). This is not the case for the afterglow evolution in a wind- like external 
medium with n oc (Chevalier &; Li 1999). Thus, a detection of the proton component in the 
GeV band lightcurve would provide a diagnostic for an approximately constant external medium. 
A non-detection, however, would not necessarily be an argument against the constant medium, 
since the phase space region for the proton component detection is small. 



6.2. Detectability of the electron IC component 

To explore the detectability of the IC component, we choose the typical parameters for regions 
II and ir to be eb = W^^eB,-2, (-e = 0.5ee,.5 and n = 1, similarly to Sari & Esin (2001). Following 
the same procedure to derive (p^, one can derive the critical time tjc when the IC component 
overtakes the synchrotron component at a typical frequency v < i'u,e- 

In the X-ray band, the overtaking time of the IC component usually occurs in the slow-cooling 
regime, which we will assume in the following discussions (see below). The crossing point between 
the synchrotron spectral component and the IC spectral component, z^g^ (Sari & Esin 2001), 
could in principle be either above or below f^'^g- If > ^Ine^ 6q.(|^) could be used directly, and 
the overtaking condition is {16/3)aT,eCenr'ym,e^^ {i^/i^c,e)^^'^ > 1- The complication in comparing 
the regions I and I' is that the Compton cooling factor (1 + Yg) in the expression of Uc,e (eq.[P) 
can no longer be neglected. For {rjee/ (-bY^"^ ^ 1, one has (1 + Yg) ~ (r/eg/e^)^''^, where in the 
slow cooling phase r/ = (7m,e/7c,e)^~^- We have derived the overtaking time for this case, where 
tic oc iy[4(4-p)/(3p2-_23p+36)] (^f^ gj^^.^ ^ ^jsin 2001). For reasonable values of p (e.g. 2.2-2.4), 
the quantity (3p^ — 23p + 36) is close to zero, causing a very sharp, and probably unphysical, 
dependence of tic on all parameters. However, the value of tic derived in this manner is not 
generally useful for our purposes here. The reason is that, to ensure fg^ > t'^g, the eg and es 
should be close to the values near the boundary between the regions II and III (eq.[^), so that 
both values are comparable. In such a regime, the approximation of (1 + yj=) ~ (r/eg/e^)^^^ no 
longer holds, and one cannot get a simple analytic expression for tic. 
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More generally, the overtaking occurs when i/g < in situations where eg e^, as is 
the case for the typical values adopted above (and in Sari &: Esin 2001). Noticing that in this 
regime Fl^,{u) = Fl^^^Ju/uj^^^y/^ = Fl%,,^,{u/uj^^,)-(P-^y'{uy^J^P-^)/^, using eq.® the 

overtaking condition can be written as (16/ 3) ax, eCenni^i/^ {vjv^^^fl'^ (i^/i^m e)^^^~^^^^ > 1- This 
finally leads to 

tic = 3.4 days e^%e^iX^(}/^e^^-V-^\X + z)-^-^^v~,l^^ (23) 

for p = 2.2 and v = lO^^Hz (X-ray band). The p-dependence in the above expression is more 
cumbersome, and is given in the Appendix. This result is in general agreement with Sari & Esin 
(2001)0, and in addition here we have explicitly presented the Cei '?52) and v dependences which 
are absent in their paper. A roughly factor of two difference on the overtaking time (~ 7.7 days in 
their case) may be caused by slightly different coefficients adopted in both works for the i'rn,e-, i^c,e, 
Fu,max,e, etc. Nonetheless, this confirms Sari & Esin's finding that in a reasonably dense medium, 
the IC component can be directly detected by Chandra a couple of days after the burst trigger. 
We note that a substantial flattening of the X-ray light curve for GRB 000926 has been detected 
by Chandra (Piro et al. 2001). Since the proton component cannot show up in the X-ray band 
under any circumstances, as we argued in §6.1, such a flattening may indicate a direct detection of 
the IC emission of the electrons ^. An alternative interpretation is advanced in Piro et al. (2001). 

Extensive efforts have been made to determine key flreball parameters such as £,€B,ee,n 
using snapshot spectral fits extending from radio to X-rays on well studied GRB afterglows (e.g. 
Galama et al. 1998; Wijers & Galama 1999; Panaitescu & Kumar 2001). A cautionary point 
which needs to be stressed about these analyses is that the spectrum which is observed need not 
be, as is generally assumed, solely due to electron synchrotron radiation, especially when using 
late time (^ a couple of days) data in the fitting. According to our results in this paper, as long 
as the CejCB phase space is in the regions I or III, the standard fitting assumption (i.e. electron 
synchrotron dominance) is safe, since there are no high energy (proton or IC) spectral components 
appearing in the X-ray band. However, in the region II of parameter space (which includes values 
of €b, Ee often derived from such fits), the analysis may not be self-consistent, since the X-ray data 
points may be due the IC component. This caution applies to even earlier snapshot spectral fits, 
if the burst happens to occur in a denser medium (notice the negative dependence of n on the 
-(2/3)(5p - 26)/(3p2 -8p- 12) index). 

The negative dependence on z/ig of tic (the explicit index is — (2/3)(3p+2)(p— 4)/(3p^ — 8p— 12), 
see also Sari &: Esin 2001) indicates that for energy bands above X-rays, the overtaking time is 



^The numerical indices in the equation (p^ ) are in agreement with Sari & Esin's resuhs for p = 2.2, except the index 
for (1 + z), where we have included an additional factor from the frequency redshift correction, i.e., -0.36=0.32-0.68. 

^After this paper was submitted, we noticed that Harrison et al (2001) have performed a detailed fit to the 
snapshot spectra of the afterglow data of this burst at 2 day and 10 day. The X-ray emission data are consistent 
with an IC emission component. The best-fit parameters they derived lie in our Regime II. This agrees with the 
discussions presented here. 
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much earliei^. In fact, in the GeV band, the IC component dominates almost throughout the 
entire afterglow phase. For the typical parameter set corresponding to the IC dominated region 
II, Figure 2b shows the time evolution of the snapshot spectra. For completeness, we present also 
in Figure 2a the snapshot spectra for the typical parameter set in the proton dominated region I, 
while Figure 2c shows the snapshot spectra for parameters in the electron synchrotron dominated 
region III. We can see that in this latter region, neither of the two high energy components is 
prominent below ~ GeV energies (although at ^ GeV there is a lower level IC component). Two 
explanations ought to be made about our code. First, to avoid adding up an unphysical component 
in the electron synchrotron self-absorption band, we have arbitrarily defined the self- absorption 
cut-off in the proton component. The self-absorption segment in the electron IC component is still 
plotted with a slope 2 rather than 1 (Sari & Esin 2001) for the convenience of code developing, 
which does not influence the final broad-band spectrum. Second, at the cut-off frequency of each 
component, we have adopted a sharp cut-off while a more realistic cut-off should be exponential. 
The same applies for the sharp jumps in the lightcurves presented in Figs. 3 and 4. 

In Figure 3, we present the X-ray lightcurves for the typical bursts in the three different 
regimes. While the regime I and III afterglows show a monotonous decay in this band, we show 
that the regime II afterglows can show interesting bump features, due to the dominance of the IC 
component at a later time, in qualitative agreement with Panaitescu & Kumar (2000) and Sari &: 
Esin (2001). 

6.3. GeV afterglows and GRB 940217 

The extended GeV emission 1.5 hours after the trigger of GRB 940217 detected by EGRET 
(Hurley et al. 1994) indicates that a high energy spectral component can extend into the GeV 
band for a long period of time, at least in some bursts. In principle, this could be either due to 
the proton synchrotron emission in the regime I, or due to the electron IC emission in the regime 
II, or even due to the electron synchrotron emission in the regime III. We will show below that the 
regime II IC-dominated origin is the more plausible explanation. 

The peak energy flux expected from various scenarios can be estimated straightforwardly. For 
the proton synchrotron in regime I, we have 

zvF,,p(GeV) = i^F,,^,,,p(i//i/^,p)-(P-i)/2 _ 14 ^ lo'^^ergs s'^ cm'^ 

xe^^-^^Ci'-^^eS^^+s'^/'n^/^fg+^^^/^I^^f (1 + z)(9-P)/4 ^-^(''-i^/Vg-^)/^ (24) 

(see also Fig. 2). This is more than one order of magnitude below the calculated level of Bottcher 
& Dermer (1998) as well as the analytical estimate of Totani (1998). The main discrepancy with 



At higher energy bands, i/ic < i^m.e may be no longer satisfied, and the overtaking time may occur in the 
fast-coohng regime. Nonetheless, the negative dependence on i/ of tic generally holds. 
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both of these results is due to their having adopted 7m,p = T) rather than the more accurate lower 
value we adopted in equation (jl]), and also due to the fact that we calculate the flux coefficient 
using p = 2.2 while they use p = 2, which further enhances the discrepancy between 7m,p and F, 
and which gives a milder i^-dependence. To test this, we have substituted 7m,p = T and p = 2 
in our code, and this reproduces Bottcher & Dermer's results. We conclude that the correction 
introduced by using here the more realistic 7m,p (eq. is essential, and that the previous rough 
estimates using 7m.p = T can considerably overestimate the proton synchrotron flux level. Another 
feature, which can be seen from equation (p^), is the negative temporal decay of the flux level with 
the index —3{p — l)/4 ~ —0.9, which indicates that even if the the proton synchrotron emission 
flux level is detectable in the early afterglow phase, it will drop with time as time goes by. For 
the regime III, electron-synchrotron-dominated case, the trend is similar, with a steeper temporal 
index -{3p - 2)/4 1.15 (Fig.4). 

In the IC-dominated regime II, contrary to the proton component, the IC component itself 
has a bump peaking at (for slow-cooling) or z^^^ (for fast-cooling) in the F^, plot. The peak 
win sweep the band u = lO^^Hz 1^23 at t = t^ax ~ 0.3hr efl^G^^'^e%^_^Sl!2^n-^l^{l + zfl^u'^'^ 
in the slow-cooling regime. The temporal index before the flux reaches its peak is 1 (slow-cooling) 
or (8 — 3p)/3(4 — p) ~ 0.26 (fast cooling), and is (11 — 9p)/8 ~ —1.1 (slow cooling, which is 
usually the case) after the flux has passed its peak. The peak energy flux can be estimated as 
^^l^max e(^ = tmax) ~ 3.0 X 10~^''ergs s~^ cm~^ for eB ~ 0.01. This temporal evolution is mild, 
which allows a substantial GeV emission component lasting hours after the GRB trigger. 

The EGRET flux sensitivity above 100 MeV is ~ 10^''ph s^^ cm^^ for point-source 
observations over a period of two weeks in directions away from the Galactic Plane. Correcting for 
an average effective on-source observing fraction of 45% (D. J. Thompson, private communication), 
the fluence threshold may be estimated as ~ lO^'^ • T{t/TY/^ ~ 5 x IQ'^t^/^ ergs cm at an 
average energy 400 MeV, where T = 14 • 86400 • 45%, and t is integration time in seconds. This 
fluence sensitivity may be extrapolated down to integration times such that at least, say, 5 photons 
are collected. For even shorter integration times, the flux sensitivity may be defined, e.g. by the 
criterion that at least 5 photons are collected, which is ~ b/{A^fft) ph s^^ cm^^. For EGRET, 
Af>ff ~ 1500cm^, giving an estimated fluence threshold in the low integration time regime of 
~ 5/1500 • 400MeV ~ 2.1 x lO'^'ergs cm^^ ^s shown in Figure 4, for £^2 ~ 1, the EGRET ^uence 
threshold level is not reached by regime I or III bursts, but this level is attainable for a regime II 
burst located at closer distances (e.g. z = 0.1), or in a higher density environment. Thus the late 
(~ hour) GRB 940217 afterglow was most likely dominated by the electron IC emission from a 
nearby or dense-medium regime II burst. This agrees with Meszaros &: Rees (1994), and Dermer 
et al. (2000b) drew a similar conclusion by detailed simulations using a specific set of parameters, 
i.e., ee = 0.5, < 10""^, which lies in our regime II. The Gamma-ray Large Area Space Telescope 
(GLAST) currently under construction will have a flux sensitivity of ~ 1.6 x 10~^^ergs s~^ cm~^ 
for long term observations (Gehrels & Michelson 1999), roughly 40 times more sensitive than 
EGRET in the point-mode. In the low flux regime, given the effective area of ~ 8000cm2 (Gehrels 
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k, Michelson 1999), it is only about a factor of 5 more sensitive than EGRET. The fluence 

threshold for GLAST is roue; hly ~ 1.2 X 10-9^1/2 

ergs cm for long integration time regime, and 
~ 4.0 X lO^^ergs cm^^ for short integration time regime (again assuming that at least 5 photons 
are collected). This will make most regime II burst afterglows detectable at a typical cosmological 
distance and in a moderate density medium (Fig.4). 

In Figure 4 we show the GeV lightcurves for bursts typical of the three different parameter 
regimes and £^2 = 1- For comparison with future observations, we have integrated over the 400 
MeV - 200 GeV band to get the total energy fluence GLAST can collect during a certain time 
duration t. The sensitivity threshold of EGRET and GLAST are indicated. One sees that regime 
II burst afterglows would be generally detectable by GL^ST within hours after the burst trigger. 
Bursts in regime I and III are generally non-detectable by GLAST, even for burst at z = 0.1. 
Increasing total energy budget {£52) or ambient density (n) can increase the detectability of these 
bursts. Since for transient events, the key factor of sensitivity is the effective collecting area, some 
future ground-based larger area GeV telescopes, such as the 5@5 - 5 GeV energy threshold array 
of imaging atmospheric Cherenkov telescopes at 5 km altitude (Aharonian et al. 2001), may have 
better chance to detect the regime I and III bursts (at the prompt phase), and of course, to collect 
more photons from regime II bursts. A criterion to differentiate between nearby regime I and III 
bursts is that the lightcurve for the regime I burst is flatter. In any case, we conclude that an 
extended GeV afterglow is a diagnostic of a regime II (IC-dominated) burst. 

Recently, an energy flux upper limit J{E > 760GeV) < 9.4 x 10~^^ergs cm~^ s~^ in GRB 
010222 was obtained in a 4-hour measurement with the stereoscopic HEGRA Cherenkov telescope 
system, 19 hours after the burst trigger (Goetting & Horns 2001). Given the cosmological distance 
of 2; > 1.477 (Jha et al. 2001), this is consistent with our model prediction in this paper (see 
Fig.2), even for the most favorable regime II case (Fig.2b). 

6.4. Prospects from broadband observations in the Swift- GLAST era 

The GLAST mission will be launched in 2005, with a sensitivity range in the 20 MeV-300 
GeV, complemented by the Glast Burst Monitor {BGM) whose energy range extends from a few 
keV to 30 GeV. Another broad band GRB mission. Swift, will be launched in 2003, and will 
be sensitive in the optical. X-ray and 7-rays up to <, 140 keV. At the same time, ground-based 
experiments such as Milagro, HESS, Veritas, MAGIG and Gangaroo-III may provide >,0.5 TeV, or 
in some cases ^30 GeV data or upper limits. In the Swift- GLAST era, simultaneous broad-band 
observations at the very earliest stages of the GRB afterglows will become possible, which will 
bring invaluable information about GRB shock physics and the central engine. Here we note 
several interesting issues which can be addressed in the Swift- GLAST era: 

1. Sari & Esin (2001) pointed out that due to the IC cooling, there are two possible solutions 
of the unknown fireball and shock parameters for a same set of observables from the low-energy 
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afterglow fits. For high density mediums, these two scenarios may be distinguished from the 
late-time X-ray observations. But for low density mediums, the two solutions are degenerate and 
indistinguishable with the present data. However, we note that the two sets of solutions lie within 
the regime III and regime II, respectively. This provides a natural way to distinguish between 
the two scenarios by using the GeV afterglow data. If an extended GeV afterglow is detected by 
GLAST, then the parameter space should be in regime II where the IC component dominates in 
the high energy band. In this case eqs.(4.17) to (4.20) of Sari & Esin (2001) will apply. Otherwise, 
the parameter space should be in regime III, where no prominent emission component shows up 
in the GeV band. This is the case that Sari & Esin's eqs.(4.13) to (4.16) may apply. 

2. It has been proposed by Waxman (1995) and Victri (1995) that GRBs are likely sites to 
produce ultra-high-energy cosmic rays (UHECRs), and that the lO^'^ eV excess UHECRs detected 
are of the GRB origin. This hypothesis is subject to debate (cf. Stecker 2000; Mannheim 2000; 
Scully & Stecker 2001). We note that future GeV observations with GLAST may be able to 
add important criteria on the debate, at least for the external shock scenario. To accelerate 
protons to ultra- high energies around 10^° eV, sb must be close to unity (Waxman 1995; Vietri 
1995; Rachen & Mcszaros 1998, but see Dermer 2001). If substantial extended GeV afterglows 
are common among GRBs, and if the redshift measurements from the low frequency afterglow 
observations indicate that the GRBs are around z ~ 1, this will impose severe constraints on the 
UHECR acceleration theory by the external shock scenario, since the regime II generally favors a 
small SB- (However, our calculations do not apply to a possible UHECR acceleration in internal 
shocks). On the other hand, if long duration GcV afterglows are not common, and for a few 
nearby bursts a GeV (prompt) lightcurve hardening is detected by GLAST or some other fTitiire 
telescopes, this could be attributable to proton synchrotron emission, in which case both €p and 
cb are close to unity. This would provide support to the theory of UHECR origin in GRBs. Since 
these components are not masked by the electron synchrotron and the IC components, this would 
also hint at a small eg (e.g. due to a weak coupling between electrons and protons, as argued by 
Totani 1998; 2000). 

3. Although most present snapshot spectral fits assume that the unknown shock parameters 
(eg, Ep, cb) are constant with time, there is no a priori reason for this simplest assumption. 
Evolutions of one or more of these parameters are in principle possible (e.g. Dermer et al. 
2000b). Present observations are too crude to explore these possibilities, but future broadband 
observations will provide the opportunity to explore this important issue. In our Figure 1, we 
present generically the parameter space of various regimes as well as their dependences on some 
other parameters including the observation time. In the whole discussions in §6, we have assumed 
the constant equipartition parameters. If there exist substantial evolutions of these parameters, 
the blastwave parameters will change in the eg — es space, and may, in some cases, switch the 
regime they belong to as time goes by. This will bring some additional interesting signatures in 
the lightcurves in various bands, and further more detailed studies may provide diagnostics on the 
possible evolutional effects. 
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7. Summary 

We have studied GRB afterglow snapshot spectra and hghtcurves over a broader band 
than usual by including the canonical electron synchrotron emission as well as two other high 
energy spectral components, i.e., the proton synchrotron emission component and the electron 
synchrotron sclf-invcrse-Compton emission component. We have in particular concentrated upon 
the X-ray to GeV-TeV ranges, including the effect of attenuation by 77 pair formation. This 
investigation has the advantage, relative to prior ones, of bringing together in a single coherent 
treatment the effects of these various high energy mechanisms, which hitherto had been mostly 
treated singly or in twos, within the context of a specific GRB afterglow dynamical model. This is 
carried out over a wider range of parameter phase space than hitherto, to allow a global view on 
the relative importance of the various spectral components. 

For the frequency range below the electron's synchrotron cut-off, v < Vu^e, there is a 
competition between the electron synchrotron component on the one hand, and the proton 
synchrotron component or the electron IC component on the other, which can affect the higher 
energy bands including X-rays or above. This competition divides the e^e^ phase space into 
three regimes (Fig.l). We have explored the range of validity of these regimes, and discussed 
the conditions for which these high energy spectral components would show up in various bands, 
especially in the GeV and the X-ray band. The conclusion is that the IC component is likely 
to be important in a relatively large region of parameter space, while the conditions for which 
the proton synchrotron component is important involve a small, but non-negligible, region of 
parameter space. One interesting consequence is that there is a substantial region (regime III) in 
which neither of the two high energy components are important. Above the electron synchrotron 
cut-off, the competition is between the electron IC component and the hadron-related photo-meson 
decay components, which we treated as a reduced extension of the proton synchrotron component. 
Again, the phase space region where the latter effects arc important in the afterglow is small. We 
also find that for the external shock and the afterglow phase the 7 — 7 absorption is not important 
below the TeV range. 

A general conclusion is that the most likely origin for an extended high energy afterglow 
component at GeV energies is from the electron IC component. Not only is the phase space 
region where the IC component dominates (II and IF) much larger than that where the proton 
component dominates (I and I', see Fig. 1), but its intensity is also much higher than that of the 
hadron components, and the time scale during which an appreciable flux level is maintained in the 
GeV band is much longer than for the hadronic components. In the parameter regime favorable 
for the IC emission, this component is observable at and above the X-ray band. In the X-ray 
band, it will lead to a flattening of the light curve at late times, as long as the medium density is 
not too low (see also Sari &: Esin 2001). Above the X-ray band, the time after which IC emission 
becomes dominant appears earlier, and the IC component dominates the GeV-band emission 
almost from the onset of the afterglow phase. In general, a high external density medium favors 
the detectability of the IC component. Such an IC component is likely to have been responsible 
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for the GeV photons detected from GRB 940217 with EGRET, and similar events should in the 
future be detectable by GLAST. 

The proton synchrotron component, as well as the hadron-related photo-meson electromagnetic 
components in the afterglow radiation are likely to be, in most cases, less important than previous 
approximate estimates indicated (Vietri 1997; Bottcher & Dermer 1998; Totani 1998; Totani 
2000). The ee,eB phase space region where this component overcomes the electron synchrotron 
and the IC components is very small, unless the medium density is high. Even in the most 
favorable region of the phase space, the proton component is not expected to show up in the X-ray 
band, but it may overcome the electron synchrotron component in the GeV band shortly after the 
onset of the afterglow phase, if the medium density is moderately high. However, the flux level 
is low and drops with time, and is generally not detectable by GLAST unless for some extreme 
bright bursts at close enough distances. We note also that the value £52 = 1 used for the estimates 
in the figures is per steradian, and narrow jets could in principle give values >,10^. The detection 



of a proton component would be of extremely high interest (§6.1) for the GRB and shock physics. 
Its detection, if successful, would imply, depending on its strength, fireball and shock parameters 
which are more extreme than currently commonly assumed. In particular, it would provide a 
diagnostic for a high Cp, and/or a low eg. 

We point out a simple way to break the current parameter space degeneracy which currently, 
from low frequency observations alone (X-rays and below), prevent the unambiguous determination 
of the unknown fireball and shock parameters, through the use of snapshot spectral measurements 
extending into the GeV range. We also suggest a way to test the hypothesis of a GRB origin for 
UHECR using combined Swift and GLAST data. Such observations may also provide diagnostics 
for the presence of a quasi-homogeneous versus a wind-like inhomogeneous external medium. 



Finally, we note that the condition for the spectral IC component to be prominent (eq.[|lf 
and Fig.l) usually covers the parameter regime in which IC cooling is important (r/eg/e^ > 1, Sari 
& Esin 2001). Thus, the current snapshot spectral fits ought to be made with caution for the 
X-ray data points when the IC cooling is important, especially for the data at later times, and for 
cases which may involve a high external medium density. 



We are grateful to C. D. Dermer, D. J. Thompson and M. Bottcher for informative 
correspondence, to the referee for valuable comments, and to NASA NAG5-9192 and NAG5-9193 
for support. 
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A. Explicit p-dependences of equations (18) and (23) 

In this appendix, we explicitly present the p-dependent indices in equations (18) and (23). 
For equation (18): 

8(P-1) 4(p-l) 8(p-2) 4(p-2) 4,3_ . 



"-4P ii_2p 4p-7 g_4p 2(p-2) 

^^^^2(16-3p)^-5(T6^^^2{16-3p)^-^^^^2Pe^^^^l6-3p _ 



For equation (23): 



3p''-10p+4 

tic = 3.4 days e Jf^^ ^H^^^ ^3i,,^-s,-i2) 

ip^-ip+l 2(5p-26) _ 3p2-16p+4 2(3p+2)(p-4) 

^^|^3p^-8p-12)^ 3(3p2_8p-12)(l + ^) 3(3p^ -8p-12) i,^g3{3p^ -8p-12) _ ^^2) 
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Fig. 1. — Regions in the €e,eB parameter space where the various radiation mechanisms dominate 
at selected frequencies. A) For photon energies ly < i'u,e (the synchrotron frequency for electrons 
at the upper end of their energy distribution) the solid lines 1 and 2 divide the space into three 
regimes. Regime I is where the proton synchrotron component overcomes the electron synchrotron 
component; regime II is where the electron IC component overcomes the electron synchrotron 
component; and regime III is where the electron synchrotron component dominates the other two. 
B) For a higher energy band with v > i^„^e, the space is divided into two regimes I and I' by a dashed 
line 3, or a dashed line 4 (4'), depending on the subcase. Regime I' is where the proton-related 
components overcome the electron IC component, and regime IF is where electron IC dominates 
over the proton components. For v < (the inverse Compton-boosted frequency of synchrotron 
photons radiated by electrons at the cooling break energy) the separation is given by line 3, which 
does not depend on the frequency. For v > the separation line is frequency-dependent, and 
given by line 4 (4'), which arc drawn for u ^ 10^6 Hz. Line 4 assumes the reduction factor is 
k = 0.1, while line 4' assumes k = 1 (see text). The dependences of the separation lines on t, n 
and 6 are indicated on the plot, which cause the different regimes to enlarge or shrink with these 
parameters. All lines are drawn using the following parameters: = = !> = 1; ^52 = 1, 
n = 1, a = 1, z = I, p = 2.2 and t = 1 hr. For a flatter p (close to 2), the lines 1, 2, 4 and 4' are 
flatter, and the regions I and I' will be slightly enlarged. For a steeper p, the trend is the opposite. 



-24- 




12 14 16 IB ZO Z2 

Log (Frequency) (Hz) 



12 14 16 IB 20 ZZ 24 

Log (Frequency) (Hz) 



-6 - 




8 10 12 14 16 IB 20 22 24 26 2S 



Log (Frequency) (Hz) 



Fig. 2. — Temporal evolution of the broad-band spectra of GRBs. Thick solid curves are the final 
spectra for various observer times, starting from (top) the onset of the afterglow, 1 minute, 1 hour, 
1 day to (bottom) 1 month, respectively. The sharpness of the breaks and cutoffs is an artifact of 
the analytical approximations, in reality these would be smoother transitions. For the top curve, 
contributions from the various radiation components are also plotted. Long dashed are electron 
synchrotron, short dashed are proton synchrotron, and dotted lines are electron inverse Compton 
emission. The thin solid line is the total energy flux level without 77 absorption correction, while 
the thick solid line is the energy flux level after the 7 — 7 self- absorption correction. The intergalactic 
absorption, which also becomes important around u = 10^^, is distance dependent and has not been 
included in this graph. Here all plots are calculated for standard parameters z = 1 (flat, A = 
universe), Ce = Cp = Ij = ^52 = 1, a = 1, z = 1, p = 2.2, and Fq = 300, while Cg, eb and n 
vary for the different regimes, (a) A typical regime I burst: eg = 10~^, = 0.5, n = 100 cm~^. 
(b) A typical regime II burst: eg = 0.5, €b = 0.01, n = 1 cm~^. (c) A typical regime III burst: 
Ee = 0.01, eB = 0.1, n = l cm~^. 
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Fig. 3. — X-ray uFy lightcurves {y = 10"*^^ Hz) for the tiiree types of bursts, starting from t^^ec 
and ending at the time when the bulk Lorentz factor F = 2. Curves I, 11, and III arc calculated 
for regimes where the proton synchrotron, the electron IC, and the electron synchrotron dominate, 
respectively. All the parameters adopted are the same as those in Fig.2, except for the regime 
II burst, where a slight denser medium n = 5cm~^ is adopted to show how the IC component 
flattens the light curve at later times. For the regime II, the contributions from the three different 
components are also shown with the same notations as in Fig.2. 
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Fig. 4. — GeV J^^^ vFi,dv ■ t lightcurves for the three types of bursts, starting from tdec and ending 
at the time when the bulk Lorentz factor F = 2. The energy flux has been integrated within the 
range of vi ~ 400 MeV to V2 ^ 200 GeV in order to compare with observations. The sohd curves 
I, II and III indicate bursts in regimes I, II, and III, respectively, at a typical cosmological distance 
= 1, for a flat A = universe). The three dotted unmarked curves are the same types of bursts 
located at 2; = 0.1. The other parameters adopted are the same as those in Fig.2. The sensitivity 
curves for EGRET and GLAST are also plotted. For low integration time regime, at least 5 photon 
detection is required. For the regime II parameters, the electron IC component gives an extended 
duration GeV emission, easily detectable by GLAST at z = I (solid line). For z = 0.1 (dotted 
lines), EGRET can detect a regime II burst, and this may account for the extended GeV detection 
from GRB 940217. 



